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1. Introduction 

The physics of ultracold atoms and Bose-Einstein condensates in optical lattices has 
made an enormous progress in the last decade. Due to the possibility to control all 
experimental parameters accurately over wide ranges and monitor the dynamics of the 
atoms in situ, optical lattices have become one of the most prominent model systems in 
quantum optics, solid state physics and nonlinear dynamics. Nowadays the shape of the 
optical potential can be engineered with astonishing precision, including in particular 
bichromatic optical lattices p~H3]. The manipulation of matter waves in the lattice 
is routinely accomplished by a static or time-dependent external field, either in an 
accelerated horizontal lattice [IJ[2J|1J[5] or a vertical lattice subject to gravity [&H9] also 
supported by magnetic levitation [7J[8]. A weak static field accelerates the atoms up to 
the edge of the Brillouin zone, where they are reflected leading to a periodic motion called 
Bloch oscillation. A strong field introduces decay by repeated Landau-Zener tunneling 
to higher bands. Thus optical lattices also became an important model system for the 
study of decay in open quantum systems (TJEJISIE] . 
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The dynamics of a quantum particle in a tilted periodic structure has been a subject 
of intensive theoretical investigations starting from Bloch's seminal paper on electrons 
in crystals [TO]. The basic structure of the spectrum and the dynamics was then clarified 
in a discussion by Zak and Wannier [TTHI3]: In particular, the spectrum is continuous 
with embedded resonance eigenstates - the so-called Wannier-Stark resonances. These 
eigenstates are arranged in ladders, where the climbing of the ladder is realized by a 
translation over one lattice period. Each ladder can be roughly associated with a Bloch 
band in the field free case. The equidistant spacing of energies in one ladder leads to 
a fully periodic motion - the celebrated Bloch oscillations. A review of these results as 
well as a surprisingly effective algorithm to calculate Wannier-Stark resonance states is 
given in [T3] , 

In the present paper we investigate a rather peculiar feature of tilted bichromatic 
lattices, the existence and properties of exceptional points. In an open system the 
eigenenergies become complex valued, where the imaginary part gives the decay rate 
(for a comprehensive discussion of such non-hermitian quantum systems see the recent 
textbook by Moiseyev [15]). If a system parameter is varied, as for example the external 
field strength, the eigenvalues show avoided crossings. For resonance states two types 
of level crossings scenarios exist - either the real parts of the energies anti cross and the 
imaginary parts cross, denoted as a type I crossing, or vice versa, a type II crossing [16J. 
A full coincidence of the real and imaginary parts is possible at isolated points in 
parameter space, the so-called exceptional points (EPs). The existence of these points 
has a remarkable implication on the dynamics of the system. Suppose the system 
parameters are varied adiabatically along a cyclic path around an exceptional point. In 
general, a cyclic evolution leaves the quantum state invariant up to a geometric phase, 
or Berry phase, which is of great importance both from a fundamental viewpoint as well 
as for applications [T71HE]. Cycling around an exceptional point has an even stronger 
effect: In addition to a geometric phase it exchanges the two crossing states. 

Both crossings of type I and II appear for Wannier-Stark systems where the periodic 
potential is sinusoidal (see, e.g., [33]), but exceptional points have not been detected. 
A modulated periodic potential with two additional parameters is, however, flexible 
enough to allow for such fascinating degeneracies and first results have been reported 
recently for a bichromatic potential [T§] . Such systems may be very well suited for 
experimental studies of the phenomena generated by such degeneracies for quantum 
systems, where experiments are still rare. 

In the following we will discuss the energy spectrum of a quantum particle in a 
tilted bichromatic optical lattice in dependence of the system parameters. A pronounced 
feature of such a Wannier-Stark system is the enhancement of the decay rate by resonant 
tunneling (RET) when the energies of two different Wannier-Stark ladders cross. Most 
interestingly, this crossing can be of both types, with an anti-crossing of either the energy 
or the decay rate, and even give rise to an exceptional point. The properties of these 
points and their embedding in the energy surfaces is discussed in detail. Exceptional 
points are most clearly identified by the divergence of the Petermann factor, which 
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measures the self-overlap of right and left eigenvectors of the same state. Thus one can 
actually use this quantity for a systematic search for exceptional points in parameter 
space. 

Finally we extend our studies to the stationary states of the nonlinear Schrodinger 
equation, which describes the dynamics of a Bose-Einstein condensate in a mean-field 
approach. The nonlinearity has a dramatic effect on the RET peaks of the decay rate, 
bending their shape and shifting their position [20]. Thus it obviously strongly influences 
the exceptional points, too. 



2. Wannier-Stark resonances 



As mentioned in the introduction, we will demonstrate the existence of exceptional 
points for a tilted periodic potential with modulated potential minima and maxima, a 
bichromatic potential. However, before we turn to the bichromatic case we first review 
some general properties of the Wannier-Stark Hamiltonian 

where the potential V(x + d) = V(x) has a single period d. It was shown [1T14T3~| |2~1~] that 
the Hamiltonian (CQ) with F > has a continuous spectrum in which discrete ladders 

£ Q ,n = E a>0 + ndF - iiy 2 (2) 

of complex valued resonances are embedded , where a = 1,2, ... is the ladder index and 
n is the site index. The decay rates T a , which are the same for all resonances in one 
ladder, are due to the finite probability for tunneling out of the lattice toward x — > — oo 
where the Stark term Fx goes to — oo (for positive fields F > 0) so that there is no 
reflection. The resonance energies (T2J) and corresponding wavefunctions thus satisfy the 
nonhermitean eigenvalue problem 

H^ a , n (x) = £ Q ,n*a,„(x) (3) 

with purely outgoing (Siegert) boundary conditions d x ^ a ^ n {x) — > —ik{x)^ a ^ n {x) for 
x — > — oo with the local wavenumber k(x) = y/2m(£ ain — V(x) — Fx)/h. At x — > oo 
the wavefunction satisfies the usual bound state condition ^^(a;)! 2 — > 0. Equivalent 
ways of avoiding reflection at x — > — oo include the use of complex absorbing potentials 
(see [15j[T9] and section H]) or complex scaling of the coordinate x [15],[22] . Here, we use 
an efficient calculation method introduced in [22],[23] based on a finite basis expansion 
in momentum space. 

The Wannier-Stark ladder (T5]) can be understood by considering the commutator 

[H, Tt] = -ldFT e (4) 

between the Wannier-Stark Hamiltonian (TjQ) and the translation operator Ti over i lattice 
sites. Equation (jlj) expresses the fact that translation by t lattice sites has the same 
effect as an energy shift IdF. Thus we obtain 

HT e V ayn {x) = TtHV^x) + [H, 7}]* a , n (a0 = (£«,„ - idF)T& ^(x) , (5) 
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which leads to the ladder (j2J) of resonance energies with the respective wavefunctions 

*a,n(z) = ^afi(x - ndF) . (6) 

The general dependence of the decay rates on the field strength is approximately 
given by T(F) oc _Fexp(— 7iAE 2 /F) where AE is the energy gap between the ground 
and first excited Bloch band. This result is obtained by Landau-Zener theory under 
the assumption that decay is mainly determined by tunneling from the ground band to 
the first excited band as successive tunneling events into higher bands, which finally 
leads to decay towards x — > — oo are fast compared to the first tunneling process 
[24"||25]. Deviations from the Landau-Zener dependence occur if a state of a lower 
ladder with energy £ Qjn is in resonance with a state of a higher ladder at a different 
site, i. e. E a ^ n = E a > tn i [SI IS]- Such a resonant coupling between two ladders leads to 
a strong increase of the decay rate of the lower ladder so that it is called resonantly 
enhanced tunneling (RET). At the same time, the decay rate in the upper ladder is 
lowered. Whenever two ladders are in resonance, the decay rate for the lower ladder 
shows a peak whereas the upper band has a dip. 

Now we turn to the bichromatic or double-periodic potential 

V(x) = y (cos(2vrx/d) + 5 cos(4nx/d + 0)) , (7) 

given by the sum of a grid with period d and an additional <i/2-periodic grid which 
creates <i-periodic potentials with modulated minima and maxima of different heights, 
depending on (f>. Throughout this paper we use scaled units such that d is equal to 2ir 
and m — H — 1. The energies are measured in units of 8Er, where Er = H 2 7r 2 / (2md 2 ) 
is the recoil energy. In the subsequent computations, we furthermore fix the potential 
strength as Vq = 1. 

The band structure of the double-periodic potential (ED) is characterized by the 
splitting of the ground band into two minibands. The Wannier-Stark ladder splits into 
two miniladders, where the width of the splitting depends on the modulation 5. In the 
limit of vanishing modulation 5 -)• the miniladders are energetically degenerate. In 
analogy to equation (J4]) one can derive the relations 

[H,T 2t ] =-2£dFT e (8) 
[H,T 2e+1 G] =-(2£+l)dFT e (9) 

(see [26l|27] for details), where H is the Hamiltonian, Ti the translation operator over 
I lattice sites and G is an operator that switches the sign of the modulation 5 in all 
following terms. Using flHJ) and ([9]) one can show that the Wannier-Stark ladder of 
eigenenergies of the unperturbed system splits into the two miniladders: 

£a, 2£ = Z a (6) + 2edF (10) 
E. a , 2e+1 = -E. a (S) + (2£+l)dF. (11) 

The energy offset E a (S) is an antisymmetric function of 5. Each Wannier-Stark ladder 
bifurcates into two which are energetically shifted with a distance 2E a (S). 
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Let us first consider the topology of the energy surfaces in an arbitrarily chosen 
part of the parameter space. Figure [1] shows the imaginary part of the eigenenergies, 
i.e. the decay rate of the two most stable states in dependence on the two parameters 4> 
and 1/F where the modulation 5 = 1 is kept fixed. 

If we vary only 1/F keeping the potential, i.e. <p, fixed, we observe the familiar 
RET spectra with eigenvalue crossings and avoided crossings. For the cut with <fi = — 1 
appearing at the front of the surfaces in figure [1] for example, we have a series of 
avoided crossings of the imaginary parts, whereas the real parts cross, i.e. a type II 
crossing scenario. Cuts at different values of <p show type I crossings. One is tempted 
to anticipate to find fully degenerate eigenvalues somewhere in the (1/F, </>)-plane. This 
will be clarified in the following section. 




Figure 1. Decay rate Ti.2 as a function of the inverse field strength 1/F and the 
phase <fi with (5=1. The configurations marked by a vertical line correspond to EPs. 

3. Exceptional points 

3.1. Overview 

Exceptional points occur in systems described by non-hermitian Hamiltonians depending 
on a set of parameters. As mentioned above, a point in parameter space at which both 
the complex eigenvalues and the eigenstates of two different eigenmodes coincide is called 
an exceptional point p~5l!28p29j . EPs must be distinguished from so-called diabolic points 
that occur if there is a coincidence of the eigenvalues but not of the eigenstates. Diabolic 
points are more familiar from hermitian systems, but in fact they are more special than 
the exceptional points and can be viewed as a coincidence of two exceptional points 
(see, e.g., [16],|30] and references therein). Bichromatic optical lattices showing diabolic 
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points in the band structure have been proposed as quantum simulators for the Dirac 
equation [3Tj . 

Before we introduce a systematic algorithm for finding the exceptional points for 
the tilted periodic potential (J7J), we first have a look at the set of solutions 8(1/ F, 5, (f>) 
shown in figure [U Exceptional points are searched by means of the criterion that the 
energy difference A£ = |£i — £ 2 | between the two lowest levels must vanish. This 
criterion is satisfied for the following points marked by vertical lines in the figure: 

1/F = 3.769, 6 = 1, = -2.991, (12) 
1/F = 6.662, 6=1, = -2.228. (13) 

One can straightforwardly verify that the two eigenfunctions involved are also degenerate 
for these parameters so that the points (TT2|) and ( fTBl are indeed exceptional points. This 
is illustrated in figure |2] which exemplarily shows the two Wannier-Stark eigenfunctions 
and ^2 for the parameters (fl2j) . One observes that the wavefunctions are localized 
in two neighboring potential minima at x = and x = 0.5. 

0.06 
0.05 
0.04 

CM 

&_0.03 
0.02 
0.01 

_5 _4 _3 _2 -1 1 2 

x/(2tt) 

Figure 2. Wannier-Stark eigenfunctions \&i and ^ in position space for the 
configuration of the EP (fT2|). The wavefunctions coincide on the scale of drawing. 




3.2. Emergence of exceptional points in tilted optical lattices 

In the following we want to address the question how the emergence of exceptional points 
in a bichromatic Wannier-Stark potential can be understood in terms of the shape of 
the potential. The definition of an EP requires the coincidence of both, the energies 
and eigenfunctions of the two states considered. 

Let us first recall the features of the single-periodic Wannier-Stark system 
introduced in section [2l where a simultaneous degeneracy of the eigenenergies and 
eigenfunctions is impossible. The Stark potential Fx leads to a localization of the 
eigenfunctions and induces a ladder structure of the eigenenergies E a>n = E aj o + 2irFn 
(see eq. ([2])). For a finite field F > the energies within a ladder differ by multiples of 
2ttF, such that a degeneracy in E can only occur between states of different ladders, 
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Figure 3. Value of the overlap integral S for a variation of the field strength F in the 
case of a single-periodic (left panel) and a double-periodic potential (right panel). The 
dashed horizontal line indicates the maximum value 5 = 1, the vertical line indicates 
the position of the maximum at 1/F = 7.5 respectively 1/F = 3.8. 



resulting in resonantly enhanced tunneling. However, the states involved are localized 
in different wells of the potential, such that they cannot be described by the same 
wavefunction. 

The situation is different in a bichromatic lattice. Within one period of the potential 
V(x), there are two potential minima, in which a Wannier-Stark resonance is localized. 
By varying one parameter, e.g. the field strength F, one can achieve a degeneracy of the 
eigenenergies. Now it is possible to fine-tune also the remaining parameters to realize a 
full coincidence as illustrated in figure HJ An exceptional point is found when the more 
stable state is destabilized and vice versa while keeping their energies degenerate. 

At the exceptional point, both the eigenenergies and the eigenstates must coincide. 
Therefore the overlap integral 

Si, 2 = |(*i|* 2 )| < 1 (14) 

of the two (normalized) eigenstates should approach unity at the EP. For parameters 
far away from it the overlap is expected to be small. As an example, figure [3] displays 
the overlap in dependence on 1/F for system parameters at some distance from an 
exceptional point (left panel), where we find a rather small peak of the overlap. The 
right panel of figure [3] displays the overlap S in dependence on F in the vicinity of the 
EP ( IT2l) . For 1/F « 3.7, which corresponds to the configuration of the EP, we observe 
a pronounced peak with S ~ 1 as the two wavefunctions are degenerate. This behavior 
can be used as a tool for detecting an EP. Note that alternatively one might consider the 
Petermann factor K a describing the '(self) over lap' between left and right eigenvectors 
that we will use in section HI 

To further investigate this phenomenon, figure H] shows both the potential and the 
wavefunction for the parameters 

1/F = 3.814, 5 = 2.251, = -3.035. (15) 
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that correspond to an EP as well as for a slightly different parameter set in the 
vicinity of the EP (see table [1] for the exact configurations). In agreement with our 
previous observations, the overlap only assumes values of S ~ 1 for a very narrow 
interval of field strengths F . The potentials shown in the two panels in figure H] are 
hardly distinguishable at first glance. Yet there is a significant difference as far as the 
corresponding pairs of Wannier-Stark eigenfunctions and in particular their overlap are 
concerned. The difference is also evident in the corresponding complex eigenenergies 
given in table [U Looking at the decay rates we observe that at the transition to the EP 



Table 1 . System parameters sets in the vicinity of the EP (|15|) . 



1/F 


S 


4> 


E x 


Ti 


E 2 


r 2 


3.000 


2.251 


-3.141 


2.561 xlO" 1 


2.961 xl0~ 2 


1.858 xlO" 1 


1.361 xlO" 1 


3.814 


2.251 


-3.035 


1.538 xlO- 1 


7.427 xl0~ 2 


1.538 xlO" 1 


7.427 xKT 2 



the ground state becomes less stable whereas the lifetime of the excited state increases. 

The emergence of the degeneracy can be related to the shape of the double-periodic 
potential. Within a period of the potential V(x) we have two minima marked by © and 
(D in the left lower panel of figure HI These minima are slightly shifted with respect 
to each other. For V® n < V® Q the ground state wavefunction is mainly localized in 
the first well (cf. left panel of figure H]). Together with well © one obtains the shape 
of a double well potential. The first excited state reveals a high probability density in 
the second well. Since the barrier in the direction x — > — oo is lower for well (D, the 
probability for tunneling processes in this direction is higher leading to a larger decay 
rate (cf. tabled]). 




-2 -1.5 -1 -0.5 0.5 1 1.5 -2 -1.5 -1 -0.5 0.5 1 1.5 

x/(2tt) x/(2tt) 



Figure 4. Configuration in the vicinity of the EP (left) and exactly at the EP ([15]) 
(right). The upper panels show the Wannier-Stark ladders for the two most stable 
resonances, where the widths of the lines indicate the decay probability. The lower 
panels show the Wannier-Stark functions V&i (blue) and ^2 (red). 
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For the field strength 1/F « 3.814 corresponding to the configuration of the EP, 
the ground and first excited resonance state are energetically fully degenerate. The first 
excited state is stabilized since tunneling towards the "more stable well" © is enhanced. 
At the same time the ground state is destabilized because the probability for tunneling 
from well © toward x —> — oo is likewise enhanced. The EP corresponds to the balanced 
situation where the two Wannier-Stark states become indistinguishable so that there 
is a degeneracy between wavefunctions belonging to different bands which cannot be 
achieved in the single-periodic Wannier-Stark system. 



3. 3. Systematic search for exceptional points 

So far all EPs have been determined in the same way: If one suspects the existence of 
an EP in some part of the configuration space, eigenvalues and eigenfunctions of the 
Hamiltonian are analyzed with respect to the conditions 

|£i-£ 2 |=0 (16) 

Sl,2=|<tfl|*2)|=l. (17) 

In the present study we numerically minimized the eigenvalues distance in ( 1T61) using 
a standard method based on the simplex algorithm (see, e. g. [32],[33]). Finally it was 
checked if the overlap condition ( fTTl) is satisfied as well. Note that a more efficient 
methods for localizing exceptional points have been developed recently [341135] . 

So far, one of the three parameters (1/F,5,<f)) has been kept fixed whereas the 
other two were varied. Now we want to search for EPs in the fully three-dimensional 
parameter space by solving the minimization problem f lT6|) . In order to converge to the 
desired EP instead of some local minimum the algorithm requires a suitable initial guess 
for the parameters (l/F,5,<p). To this end we represent the parameter space in terms 
of spherical coordinates 

A(l/F) = r sin d cos ip 

AS = r sin •& sin tp (18) 
A(f) = r cos'd. 

If the parameters r EP = (1/F EP , <5 EP , EP ) of an EP are known, one can look for 
further solutions r EP = r EP + Ar in its vicinity. Thus we search for EPs in the following 
manner: 

(i) Find an EP (e. g. with the methods described further above); 

(ii) Choose some small "distance" r, in which to look for the next EP (see below); 

(iii) Minimize the condition ( Ti~6l) using the simplex algorithm with the parametrization 
(r = const, if). 

(iv) Repeat this procedure starting with the configuration of the newly found EP. 
For the sake of completeness we note that we use 

r(A(l/F), A5, A0) = V(A(1/F))2 + (A5)2 + (A0)2 (19) 
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with A(l/F) = 1/Fi — l/F 2 , A8 = 81 — 82 and A</> = (pi — 4> 2 as the distance between 
two configurations [l/F 1 , 81, 4>i) and (l/F 2 , 8 2 , <fi 2 ). 

With this procedure the minimization problem ( lT6i) is solved in a two-dimensional 
subspace (the surface of a sphere) of the three dimensional parameter space which 
reduces the numerical effort. The finite area of that subspace also improves the 
convergence. This systematic search detects all exceptional points between the two 
ladders in the parameter region under consideration. Obviously, further exceptional 
points can exist between higher ladders, which are not detected. However, these are of 
less practical interest, as the decay is significantly stronger. 
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Figure 5. Left panel: Location of the EPs in parameter space for the two most stable 
resonances (red curve). On the plane for S = 1 the absolute value of the difference of 
the decay rates AT = \Ti — L2I is shown as a contour plot. Right panel: Curve traced 
out by the EPs projected onto the (l/F, <5)-plane. 

Starting with the EP (112jl we obtain a one-dimensional manifold formed by the EPs 
of the two most stable resonances of the double-periodic Wannier-Stark system that is 
embedded in the three-dimensional parameter space. This is displayed in figure |5j In 
addition, we have included a contour plot of the absolute difference Ar = |Ti — r 2 1 of 
the decay rates on the plane 8 = 1. The intersections of the EP curve with this plane 
are located at the zeros of ( fT6l) . For 8=1, three solutions can be detected in the region 
considered. 

One can draw the following conclusions: In the single-periodic system no 
exceptional points were found. The results for the bichromatic system reveal the 
transition between the two different potentials: In the limit 8 — > the position of 
the EP is shifted toward an infinitely strong field strength F (l/F — > 0), so that an EP 
cannot exist in the single-periodic limit. 

Finally we briefly want to discuss the shape of the curve. From figure |5] we see 
that the EP curve bends strongly for 8 ~ 2.3, so that no EPs can be found for higher 
values of the modulation 8. Yet for small field strengths F, i. e. for high values of l/F, 
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Figure 6. Cyclic variation of the parameters F and (j> for <5 = 1. The path © encloses 
an EP ((g)), the other paths do not. 

there are intersections of the EP curve with the planes (1/F,5 = const, 0). This local 
behavior repeats itself and can also be observed for higher bands. The results in the 
domain of interest support our expectation that there are no EPs in the limit F — > 0, 
i.e. without an external field. To gain further insight into the characteristic shape of 
the EP curve additional investigations are needed. 

3.4- Cyclic parameter variation 

In the following we want to study the behavior of eigenvalues and eigenvectors under 
cyclic parameter variation. As an example we consider the vicinity of the EP ( 1121) and 
choose 6 = 1 to be constant. The quantities F and are varied along the paths 

l/F = l/F EP + rsin(/3) + l/F 

= 0ep + r cos(/3) + 0o 

parametrized by an angle G [0,27r). As in subsection 13.31 r is a numerical quantity 
that can be identified as a radius in parameter space. 

Apart from paths that enclose the EP (1/F EP , Sep, 0ep) we are also interested in 
paths that do not enclose it, which can be controlled via the shifts F Q and O . Figure 
[6] shows nine selected paths, the corresponding eigenvalue trajectories in the complex 
energy plane are displayed in figure [71 

We observe two qualitatively different kinds of behavior: Any closed path in 
parameter space that does not enclose the EP leads to eigenvalue trajectories in the 
complex energy plane that consist of two separated closed curves, one for the ground 
state and one for the first excited state. Thus the eigenvalues always remain on the 
same Riemann sheet. A different result is obtained for the path © that encloses the EP. 
Instead of two separate eigenvalue curves in the complex energy plane we now find a 
single closed curve consisting of the complex eigenvalues of both eigenstates involved. 
This property is typical of exceptional points [36] and provides a useful criterion for 
proving the existence of an exceptional point in experiments (see e. g. j37J[38])- After 



Exceptional points in bichromatic Wannier-Stark systems 



12 



two full cycles the initial eigenvalue is recovered, i.e. the eigenvalues vary with a period 
of 4-7T in parameter space. 

It is also instructive to look at the eight outer subplots in figure [7] in the cyclic 
sequence 

©->-®->-®->-©->-®->-®->-®->-©->-©. (21) 

Comparing the eigenvalue trajectories for path © and path © we observe that the 
assignment of the two trajectories to the ground state (blue) and the first excited state 
(red) is interchanged. This can be interpreted as as closed path in parameter space that 
also encloses the EP and hence moves from one Riemann sheet to the other. 

Apart from the eigenvalues we are also interested in the eigenfunctions. In order 
to monitor their behavior, it is sufficient to consider only one complex number, namely 
the projection of each wavefunction on one of the basis vectors in the plane-wave basis 
chosen for the calculations. Here we choose the projection (\£ / Q ,) ma x with the maximum 
value of |^/a| 2 ; the other projections show a qualitatively similar behavior. In figureElwe 
display the phase Q (left panel) as well as the real and imaginary parts (right panel) of 
(f a ) max = |(\l/a)max| e~ l2lTn for the ground and first excited state for a cyclic parameter 
variation around the EP, i.e. when the phase /3 in parameter space varies form to 
2tt. The figures reveal the expected behavior: After one complete cycle around the EP 
the phases O of the two components are interchanged, however with a phase change 
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Figure 7. Eigenvalue trajectories in the complex energy plane £ = E — iF/2 for a 
cyclic parameter variation of F and (f> for 5 = 1 , cf. also figure [51 The energy of the 
EP is marked by ®. 
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of 7r, i.e. a different sign. In the right panel, showing the eigenvector components in 
the complex plane, one observes the resulting 87r-periodicity. Starting at the point (a) 
(ground state) and (d) (first excited state) for ir = 0, the selected components (^i)j 
reaches points (b) and (a) after one cycle, where (b) corresponds to the sign-changed 
component (^2)1 at (d). The dashed lines show the development for three further cycles. 

The periodicity of the eigenvalues and eigenf unctions can be summarized in the 
following diagrams (cf. [16] ) 




where ^1, ^2 are the Wannier-Stark eigenstates involved in the formation of the EP 
and £1, £2 the corresponding eigenenergies. 

Finally, looking again at figure we observe that in the vicinity of the parameter 
values 1/F = 4, 5 = 2.3 and <p = —3 two exceptional points approach each other. 
This region is magnified for 5 = 2.3 in figure [9], which shows the energy difference 
A£ = |£i — £ 2 | in dependence on 1/F and 0. It is instructive to have a closer look 
at the surfaces of the real and imaginary part of the eigenenergies in parameter space 
in this region. This is illustrated in figure [101 which shows the Riemann surfaces of 
the complex energy in the vicinity of these two EPs. The observed local behavior 
resembles the "double-coffee-filter" scenario typical of unfolding a diabolic point into 
two EPs [T6 |[30l[39H4T] . Note that the apparent gaps in the Riemann surfaces are an 
artifact of the finite discretization of the parameter space. 




Figure 8. Component {^ a )m&x = ('l'a)max e~ i2nQ of the ground state (blue) and the 
excited state (red! for a cyclic parameter variation around the EP (Tl2l) . Shown is the 
phase fl (left panel) as well as the real and imaginary parts (right panel). 
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-3 




l/F 

Figure 9. Difference of the complex energies | £2(1 / i 7 ^ 4>) — 4>)\ as a contour 

plot for 6 = 2.3. The marked minima correspond to the configuration of the EPs. The 
absolute values of the energy difference are color coded. Blue indicates for zero, red 
the maximum value. 

3. 5. The Petermann factor 

For a non-hermitian system, one has to distinguish between right and left Wannier-Stark 
states, i.e. the eigenstates of H for eigenvalue E a and of W for eigenvalue E^, denoted 
as \^ a ) and \^ a ), respectively. These states form a bi-orthogonal set, i.e. (^f a \^ a ') = 
for a 7^ a 1 . At an exceptional point, two eigenstates of H coincide, as well as the 
corresponding eigenstates of W , with the consequence that the scalar product {^/ a \^ a ) 
vanishes. A useful quantity in this context is the Petermann factor [4"2||4"3] 





|*a>(*a 










2 



(24) 



which provides a measure of the overlap between left and right eigenvectors. Note 
that K a is independent of the normalization of the states and diverges at the EP. The 
Petermann factor provides a convenient tool for localizing an exceptional point. 




Figure 10. Energy surfaces £1, £2 in the real- (left panel) and imaginary part (right 
panel) for 8 = 2.3. 
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Figure 11. Petermann factor K a of the ground state {K\, blue) and the first excited 
state {K 2, red) for a double-periodic Wannier-Stark system. The parameters are in 
the vicinity of the EP (fl2]l. 



Figure [IT] shows the Petermann factor of the ground state and the first excited 
state in dependence on the field strength F . The modulation 5 and phase <j) were chosen 
according to the configuration of the EP ffTiZ]) . The offset between the two curves is 
caused by the difference in the decay rates of the two respective states. At the position 
corresponding to the EP both Petermann factors show pronounced peaks. Due to 
the discretization of the field strength and the fact that the configuration (fT2|) is only 
very close to but not exactly in the EP, the numerical value of K a remains, of course, 
finite. 

4. Nonlinear crossing scenarios 

One of the most promising realization of Wannier-Stark systems, both theoretically 
and experimentally, is based on atomic Bose-Einstein condensates in tilted optical 
lattices. Near the zero temperature limit these systems can be described by a nonlinear 
Schrodinger equation (NLSE) or Gross-Pitaevskii equation (GPE) (see e. g. [HHU]) 



in a mean- field approach, where the nonlinear term g|\I/(:r)| 2 takes into account the 
interaction between the condensate particles. The interaction between the particles 
can either be repulsive (g > 0) or attractive (g < 0). In the case of Wannier-Stark 
resonances, the chemical potential /i = M — iT/2 is complex where the real part M 
describes the position of the resonance and the imaginary part T/2 accounts for the 
decay rate. 

Up to now, complex eigenvalue crossing scenarios and exceptional points in the 
context of nonlinear systems have mostly been investigated by means of simple nonlinear 
non-hermitian two-level models [4*8H50] . three- level systems [51], or to some extent 



~_tf_<¥_ 
2m dx 2 



+ V(x) + Fx + g\^(x)\ 2 ^(x) = n^(x) 



(25) 
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Figure 12. Crossing behavior of the chemical potential fi = M — lT/2 of the two most 
stable resonances for g = —0.02 (blue), (black) and 0.02 (red) in the vicinity of the 
EP (TT2]) for g = 0, where real parts (left) and imaginary parts (right) are degenerate 



at F = 0.2653. 




analytically solvable model systems like one- dimensional 5- and 5-shell potentials [52] 

or multiple-barriers [53] of a delta-comb [M] • An exception is the discovery of EPs for 
atomic gases with an attractive 1 / r- interaction reported in [38] . 

The numerical calculation of nonlinear Wannier-Stark resonances is a nontrivial task 
because most methods established for the linear system cannot be straightforwardly 
adapted since the superposition principle is no longer valid. In [53, [56] nonlinear 
Wannier-Stark resonances of the ground ladder in single-periodic and double-periodic 
systems were calculated using a complex scaling procedure developed in [57],[58]. Here 
we use a method based on complex absorbing potentials in position space [19] as it is 
the only procedure so far that has successfully been used to calculate Wannier-Stark 
resonances of excited ladders. In the following we want to briefly discuss the impact of 
the nonlinear interaction on eigenvalue crossing scenarios and exceptional points for the 
bichromatic WS-system discussed in the preceding sections. 

As a first example, figure [12] shows the real and imaginary parts of the eigenvalues 
of the chemical potential \x in dependence on the static field strength F for a WS-system 
again in the vicinity of the EP ( fl~2l . The black curves are the real and imaginary parts 
of the resonance energies of the two most stable states for the non-interacting system 
(g = 0) which cross at F = 0.2653. A moderate nonlinearity removes this degeneracy 
and we observe a type I crossing (cf. section [1]) for a repulsive interaction (g = +0.02, 
red curves) and a type II crossing for an attractive interaction {g = —0.02, blue curves) 
in agreement with the results already briefly reported in [19] . 

This indicates that the position of the EP is shifted by the nonlinearity in opposite 
directions for positive or negative values of g, assuming that the existence of the EP 
is not affected by weak interactions. Figure [13] displays the Petermann factors K\ and 
Ki for the resonance states shown in figure [12] The Petermann peaks for g = ±0.02 
are now strongly reduced and shifted to stronger fields for repulsive and weaker fields 
for attractive interactions. Moreover, the two maxima are also shifted relative to each 
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Figure 13. Pctcrmann factor of the ground state K\ and of the first excited state K-i 
for a parameter variation of F around the EP (|12[) of the linear system. The interaction 
strength was chosen as g = —0.02 (blue),0 (black) and 0.02 (red). 

other: The maxima of the excited state are more affected by the nonlinearity than the 
ground state. 

Finally we have a closer look at the range of moderately increased attractive 
interaction strengths g with —0.3 < g < 0. The corresponding results are shown in 
figure HH for the ground state. With increasing \g\ we observe that the Petermann 
peaks get more and more "blurred" and become strongly asymmetric, indicating that 
the system moves away from the true degeneracy in parameter space. 

An extension of the systematic search for EPs as described in section l3~3l to nonlinear 
systems deserves clearly further studies. In particular, it has to be checked if the 
exceptional point exists for all values of g, or if it vanishes when the nonlinearity exceeds 
a critical value. However, such studies are computationally demanding, in particular 




0.26 0.265 0.27 

F 



Figure 14. Petermann factor K 1 of the ground state for a variation of F around the 
configuration of the EP (|12|) in the linear case. The interaction strength is chosen as 
g = 0, -0.1 (blue), -0.2 (red) and -0.3 (magenta). 
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because of the convergence properties of the method for calculating nonlinear Wannier- 
Stark resonances. Thus substantial modifications of the presented numerical procedures 
are required. In this context the criterion of a vanishing inverse Petermann factor 1 / K a 
for the configuration of an EP might help. One advantage of this criterion is that it 
only requires the eigenstates of a single ladder. 

5. Conclusion 

This paper investigated the spectrum of an experimentally realizable model system for 
decay in open quantum systems, namely a quantum particle in a tilted bichromatic 
lattice. The analysis concentrated on the exceptional points (EPs) of the system, i.e. 
points in parameter space where both the complex eigenvalues and the eigenf unctions 
of two different eigenmodes coincide, which had so far mostly been considered for very 
simple toy models. An efficient method for finding EPs was presented and their location 
within the parameter space was discussed in detail as well as the properties of the 
corresponding degenerate eigenfunctions. It was demonstrated that in the limit case 
of a monochromatic lattice there can be no exceptional points for finite values of the 
tilt, because degeneracies of the eigenenergies can only occur for eigenstates of different 
ladders which cannot coincide as they localize in different lattice sites. Furthermore 
the geometric phases (Berry phases) and the eigenvalue trajectories in the complex 
plane for closed paths in parameter space were considered. The eigenvalue trajectories 
show the familiar behavior; they form a single closed curve for paths in parameter 
space that enclose an EP but two distinct closed curves otherwise. In particular, the 
two crossing states are interchanged after a full cycle in parameter space enclosing an 
EP. Finally the case of an interacting Bose-Einstein condensate in a tilted bichromatic 
optical lattice was discussed within a mean-field approach by solving the corresponding 
nonlinear Schrodinger equation or Gross-Pitaevskii equation. It was found that the 
type of crossing, i. e. if the real parts of the eigenenergy anti-cross while the imaginary 
parts (decay rates) cross (Type I) or vice versa (Type II), depends on the sign of the 
interaction between the particles (i.e. whether it is attractive or repulsive). These 
crossing scenarios, as well as an analysis of the Petermann factor, which measures the 
overlap between left and right eigenstates and shows a pronounced peak in the vicinity 
of an EP, indicates that the EPs are shifted by the nonlinear interaction in a direction 
opposite to its sign. 

Acknowledgments 

We would like to thank Sandro Wimberger for fruitful discussions. Financial 
support by the Deutsche Forschungsgemeinschaft (DFG) via the Graduiertenkolleg 
792 "Nichtlineare Optik und Ultrakurzzeitphysik" is gratefully acknowledged. K.R. 
acknowledges support from a scholarship of the Universit Libre de Bruxelles. 



Exceptional points in bichromatic Wannier-Stark systems 



19 



References 

[1] T. Salger, C. Geckeler, S. Kling, and M. Weitz, Phys. Rev. Lett. 99 (2007) 190405 
[2] G. Ritt, C. Geckeler, T. Salger, G. Cennini, and M. Weitz, Phys. Rev. A 74 (2006) 063622 
[3] S. Foiling, S. Trotzky, P. Cheinet, M. Feldand R. Saers, A. Widera, T. Miiller, and I. Bloch, Nature 
448 (2007) 1029 

[4] M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, Phys. Rev. Lett. 76 (1996) 4508 
[5] C. Sias, A. Zenesini, H. Lignier, S. Wimberger D. Ciampini, O. Morsch, and E. Arimondo, Phys. 

Rev. Lett. 98 (2007) 120403 
[6] B. P. Anderson and M. A. Kasevich, Science 282 (1998) 1686 

[7] M. Gustavsson, E. Haller, M. J. Mark, J. G. Danzl, G. Rojas-Kopeinig, and H.-C. Nagerl, Phys. 

Rev. Lett. 100 (2008) 080404 
[8] M. Gustavsson, E. Haller, M. J. Mark, J. G. Danzl, R. Hart, A. J. Daley, and H.-C. Nagerl, New 

J. Phys. 12 (2010) 065029 
[9] Q. Beaufils, G. Tackmann, X. Wang, B. Pelle, S. Pelisson, P. Wolf, and F. Pereira dos Santos, 
Phys. Rev. Lett. 106 (2011) 213002 
[10] F. Bloch, Z. Phys 52 (1928) 555 
[11] J. Zak, Phys. Rev. Lett. 20 (1968) 1477 
[12] J. Zak, Phys. Rev. 181 (1969) 1366 
[13] G. H. Wannier, Phys. Rev. 181 (1969) 1364 

[14] M. Gluck, A. R. Kolovsky, and H. J. Korsch, Phys. Rep. 366 (2002) 103 

[15] N. Moiseyev, Non-Hermitian Quantum Mechanics, Cambridge University Press, Cambridge, 2011 
[16] F. Keck, H. J. Korsch, and S. Mossmann, J. Phys. A 36 (2003) 2125 

[17] A. Shapere and F. Wilczek, Geometric Phases in Physics volume 5 of Advanced Series in 

Mathematical Physics, World Scientific, Singapore, New Jersey, London, Hong Kong, 1989 
[18] E. Sjoqvist, Physics 1 (2008) 35 

[19] K. Rapcdius, C. Elsen, D. Witthaut, S. Wimberger, and H. J. Korsch, Phys. Rev. A 82 (2010) 
063601 

[20] K. Rapedius and H. J. Korsch, Phys. Rev. A 77 (2008) 063610 

[21] J. E. Avron, J. Zak, A. Grossmann, and L. Gunther, J. Math. Phys. 18 (1977) 918 

[22] M. Gluck, A. R. Kolovsky, H. J. Korsch, and N. Moiseyev, Eur. Phys. J. D 4 (1998) 239 

[23] M. Gluck, A. R. Kolovsky, and H. J. Korsch, J. Phys. A 32 (1999) L49 

[24] M. Holthaus, J. Opt. B 2 (2000) 589 

[25] M. Gluck, A. R. Kolovsky, and H. J. Korsch, J. Opt. B: Quantum Semiclassic. Opt. 2 (2000) 694 

[26] B. M. Brcid, D. Witthaut, and H. J. Korsch, New J. Phys. 8 (2006) 110 

[27] B. M. Brcid, D. Witthaut, and H. J. Korsch, New J. Phys. 9 (2007) 62 

[28] T. Kato, Perturbation theory for linear operators, Springer Verlag, Berlin, 1966 

[29] W. D. Hciss, Phys. Rev. E 61 (2000) 929 

[30] M. V. Berry and M. R. Dennis, Proc. Roy. Soc. Lond. A 459 (2003) 1261 
[31] D. Witthaut, T. Salger, S. Kling, C. Grossert, and M. Weitz, arXiv:1102.4047 (2011) 
[32] D. J. Wilde and C. S. Beightler, Foundations of Optimization, Prentice-Hall, Inc., Englcwood 
Cliffs, N.J., 1967 

[33] J. C. Lagarias, J. A. Reeds, M. H. Wright, and P. E. Wright, SIAM Journal on Optimization 9 
(1998) 112 

[34] R. Lefebvre and N. Moiseyev, J. Phys. B 43 (2010) 095401 
[35] R. Uzdin and R. Lefebvre, J. Phys. B 43 (2010) 235004 
[36] M. V. Berry, Proc. R. Soc. Lond A 392 (1984) 45 

[37] H. Cartarius, J. Main, and G. Wunner, Phys. Rev. Lett. 99 (2007) 173003 
[38] H. Cartarius, J. Main, and G. Wunner, Phys. Rev. A 77 (2008) 013618 
[39] A. P. Scyranian, O. N. Kirillov, and A. A. Mailybaev, J. Phys. A 38 (2005) 1723 
[40] O. N. Kirillov, A. A. Mailybaev, and A. P. Scyranian, J. Phys. A 38 (2005) 5531 



Exceptional points in bichromatic Wannier-Stark systems 



20 



[41] O. N. Kirillov, J. Phys.: Conf. Ser. 181 (2009) 012023 

[42] K. Petcrmann, IEEE J. Quantum Electron QE-15 (1979) 566 

[43] M. V. Berry, J. Mod. Opt. 50 (2003) 63 

[44] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation, Oxford University Press, Oxford, 2003 

[45] A. S. Parkins and D. F. Walls, Phys. Rep. 303 (1998) 1 

[46] F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71 (1999) 463 

[47] A. J. Leggctt, Rev. Mod. Phys. 73 (2001) 307 

[48] E. M. Graefe and H. J. Korsch, Czech. J. Phys. 56 (2006) 1007 

[49] D. Witthaut, E. M. Graefe, and H. J. Korsch, Phys. Rev. A 73 (2006) 063609 

[50] E. M. Graefe, U. Giinthcr, H. J. Korsch, and A. E. Niederle, J. Phys. A 41 (2008) 255206 

[51] E. M. Graefe, H. J. Korsch, and D. Witthaut, Phys. Rev. A 73 (2006) 013617 

[52] D. Witthaut, S. Mossmann, and H. J. Korsch, J. Phys. A 38 (2005) 1777 

[53] K. Rapedius and H. J. Korsch, J. Phys. B 42 (2009) 044005 

[54] D. Witthaut, K. Rapedius, and H. J. Korsch, J. Nonlin. Math. Phys. 16 (2009) 207 

[55] S. Wimberger, P. Schlagheck, and R. Mannella, J. Phys. B 39 (2006) 729 

[56] D. Witthaut, E. M. Graefe, S. Wimberger, and H. J. Korsch, Phys. Rev. A 75 (2007) 013617 

[57] N. Moiseyev and L. S. Cederbaum, Phys. Rev. A 72 (2005) 033605 

[58] P. Schlagheck and T. Paul, Phys. Rev. A 73 (2006) 023619 



